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. ABSTRACT 

o . 

We present a new set of spectroscopic and photometric data extending the obser- 
vations of SN 1997D to over 400 days after the explosion. These observations confirm 
the peculiar properties of SN 1997D, such as the very low abundance of 56 Co (0.002 
r- 1 ! M©) and the low expansion velocity of the ejecta (~ 1000 km s _1 ). We discuss the 

O implications of these observations for the character of the progenitor and the nature 

of the remnant showing that a Crab-like pulsar or an accreting neutron star formed 
in the explosion of a low mass progenitor should already have produced a detectable 
luminosity at this epoch, in contrast with photometric data. On the other hand, the 
explosion of a high mass progenitor with the formation of a black hole is consistent 
with the available observations. The consequences of this conclusion regarding the na- 
ture of the explosion and the prospects of directly identifying the black hole are also 
addressed. 
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1 INTRODUCTION 

SN 1997D is the least luminous and least energetic type II 
supernova discovered to date. An earlier analysis (Turatto 
et al. 1998, hereafter Paperl) showed that the SN peak mag- 
nitude was fainter than Mb= —14.65, had a very red color, 
indicating that the photospheric temperature at discovery 
was only T e // = 6400° K, and a very slow expansion veloci- 
ties of the order of 1000 kms -1 . The luminosity in the late 
light curve indicated that the ejected 56 Ni mass was only 
O.OO2M0, at least one order of magnitude smaller than the 
estimates for normal type II supernova (SNII). 

A similar small amount of 56 Ni was suggested for an - 
other Type II supernova, SN 1994W flSollerman et al,1998[). 



Based on observations collected at ESO-La Silla (Chile) 



However Sollerman et al. (1998) argues that the light curve 
of SN 1994W is contaminated in different phases by the 
probable contribution of CSM-ejecta interaction and by dust 
formation in the ejecta which makes this case less compelling 
than that of SN 1997D. 

The unusual properties make SN 1997D an intriguing 
target for further study, especially with regards to the char- 
acter of the progenitor and the nature of the compact rem- 
nant left after the explosion. In the following we will discuss 
two extreme alternatives: a high mass (25-40M©) star in 
which the low mass of 56 Ni is th e result of the fallback o f 
material onto the collapsed core ( [Woosley fc Weaver 1995| ), 
or a low mass (8-IOM0) star in which hardly any 00 Ni 
is produced to begin with (Nomoto et al. 1982; Nomoto 
1984, 1987). In Paperl, based on modeling of the light curve 
and spectrum, we favored a scenario in which a high mass 
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(26M©) progenitor exploded about 50 days before the dis- 
covery with an explosion energy of only 4 x 10 50 erg. Instead, 
Chugai & Utrobin (2000) using a hydrodynamical model and 
an analysis of the nebular spectrum suggested that a low 
mass progenitor (and a low explosion energy) is preferred. 
Identifying which of the two scenarios lies behind SN 1997D 
is of significant interest, as well as being an observational 
challenge. In particular, there would be several intriguing 
consequences if the progenitor was indeed a high mass star, 
which would indicate the existence of a novel sub-class of 
low-energy SNII. 

Crucial to our discussion is the fact that the two sce- 
narios differ qualitatively with regard to the nature of the 
compact remnant they produce. In the case of a low mass 



Table 1. Main data of SN 1997D and NGC 1536 



progenitor, the likely remnant is a neutron star ( Nfomoto et 
al. 198^), while for a high mass progenitor the rate of early 
fallback during the explosion is expected to induce the col- 
lapse of the proto-neutron star into a black hole (Colgate 
1971; Chevalier 1989; Woosley & Weaver 1995; Woosley & 
Timmes 1996; Fryer 1999). The supernova birth of a black 
hole has recently found an important indirect observational 
proof from the detection of an overabundance of cv-elements 
in the atmosphere of the companion star orbiting a prob- 
able black hole in th e X-ray binary system GRO J 1655-40 
([sraelian et al. 1999). 



A direct observation of the effect of the newly formed 
black hole on the light curve is usually impossible, owing 
to the large abundance of radioactive isotopes. Zampieri et 
al. (1998b) found that if SN 1987A had produced a black 
hole, it would emerge as late as 900 years after the explosion 
at a luminosity of ~ 10 32 erg s" 1 . In the case of SN 1997D, 
Zampieri et al. (1998a) pointed out that because of the very 
low absolute magnitude due to the low abundance of ra- 
dioactive isotopes in the ejecta, the contribution by late time 
accretion of material onto the black hole could emerge above 
the radioactive decay about three years after the explosion. 
Since this luminosity is expected to have a characteristic de- 
cay in time of t -25 / 18 (rather than the exponential decline 
typical of radioactive sources), it should leave a distinguish- 
able imprint in the late time light curve. This estimate was 
recently confirmed by a detailed investigation by Balberg 
et al. (2000), who found that the luminosity at emergence 
should be 0.5 — 3 x 10 36 erg s _1 . Such a luminosity, at the dis- 
tance of SN 1997D, is marginally detectable with HST, so, 
in principle, SN 1997D could offer a first direct observation 
of a black hole in a supernova. 

In this work, we revisit the unusual character of 
SN 1997D and present the full set of photometric and spec- 
troscopic data up to about 400 days after the explosion. In 
§ 3 we discuss the decline of the late time light curve in con- 
nection with the nature of the progenitor of SN 1997D and 
its compact remnant and show that the presence of a black 
hole is consistent with the available data. We present our 
conclusions and discuss several implications that arise from 
them in § 4. In particular, we emphasize the importance of 
extending the observations of SN 1997D to the detection 
limit of the available ground and space based instruments. 



2 THE AVAILABLE DATA SET 



Parent galaxy 
Galaxy type 

^ helio 

Distance modulus (Ho 
Galactic Ag 

SN type 

magnitude at max 
RA(SN) (2000) 
Dec(SN) (2000) 
Offset from nucleus 
Date of discovery 
Date of maximum 
h bo i at maximum 



NGC 1536 
SB(s)c pec: t 
1296 kms-'t 
75) 30.64* 
0.0* 

II peculiar 

Vmax *C 16 

04 h ll m 01.0 s * 

-56°29'56"0 * 

11" E 43" S * 

1997 Jan. 14.15 UT 

JD ~ 2450430 ± 20 (1996 Dec. 12) 

> 1.23 X 10 41 erg s" 1 



t RC3; t Tully (1988); * Barbon et al. 1999 



01: 




100 200 300 400 500 600 700 
JD-2450430 

Figure 1. B, V, R, and I light curves of SN 1997D. Three of the 
prediscovery upper limits reported in Paperl are shown. 



2.1 Photometry 

The main data relative to SN 1997D and to the parent 
galaxy are summarized in Tab. [j] (cfr. Paperl). 

All the photometric measurements available to date are 
listed in Tab. ^. Observations were obtained on 27 different 
epochs up to about 400 days after explosion. Data reduc- 
tion followed the standard procedures making use for the SN 
measurements, of a PSF fitting technique. The mean pho- 
tometric errors, estimated with artificial stars experiments, 
are about 0.05 mag when the SN was brighter than ~ 19.5 
mag, 0.10 mag for magnitudes ranging between 19.5—21 and 
0.2 mag for the faintest magnitudes. 

The B, V, R and I light curves of SN 1997D are shown 
in Fig. jjj. On the basis of the light curve modeling, in Pa- 
perl we argued that the maximum occurred around J.D. 
2450430 (±20d) and that the SN has been caught at the end 
of the plateau phase on the way to reach the radioactive tail. 
With the new data the late time light curve appears to show 
two distinct behaviors: from day 70 to about 200 the de- 
cline rate 7 is much steeper at red wavelengths than at blue 
ones (tb = 0.16 mag(lOOd)" 1 , 7v = 0.78 mag(lOOd) -1 , 
7H, = 0.80 mag(100d)~ 1 and 71 = 0.96 mag (100d) _1 ), while 
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Table 2. Photometry of SN 1997D 
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Figure 2. Bolomotric light curve of SN 1997D. The filled circle 
includes all the flux between 0.35 and 2.5/i. Open circles represent 
the BVRI flux scaled to match the filled circle at the common 
epoch. The 56 Co decay line is also reported for comparison. 



after day ~200 the slopes are similar to each other (73 = 
0.91 mag(lOOd)" 1 , 7v = 0.99 mag(lOOd) -1 , 7r = 1.00 
mag (100d) _1 , 71 = 1.18 mag (100d) _1 )and remarkably close 
to that of 5e Co (7 = 0.98 mag (lOOd)" 1 ). 

We used the multicolor photometry to estimate the 
bolometric light curve shown in Fig. ^. This has been ob- 
tained by integrating the flux in the BVRI bands, and then 
scaling up the curve to match the spectrophotometric obser- 
vation of Jan. 27-31 (filled symbol in Fig. ^ spanning the 
range 0.35 to 2.5 fj.. This means we made the assumption 
that the IR contribution to the total luminosity remains 
constant at la ter epochs, consistently with the finding of 
Schmidt ( 2000). We note that most of the flux is emitted in 
the V, R and I bands and that, owing to the slow decline, the 
contribution of the B band becomes significant only after 250 
days. Looking in more detail we find that the decline of the 
bolometric luminosity in the period 60-200 days, 0.89 ±0.02 
mag(lOOd) -1 , is marginally slower and instead, after day 
200, 7boi = 1.07 ±0.03 mag(lOOd)" 1 , is slightly higher than 
that expected if the energy is supplied by the decay of Co. 

The ratio of the bolometric luminosities of the last 
segment of the light curve with the corresponding one of 
SN 1987A is L(87A)/L(97D) = 32.5, somewhat smaller than 
the value derived early on (Paperl). Assuming for SN 1997D 
the same 7— deposition as in SN 1987A, we confirm our ear- 
lier suggestion that the mass of ejected 56 Ni is about 0.002 
M . 



2.2 Spectroscopy 

The journal of the spectroscopic observations is given in 
TabQ Figure ^| illustrates the spectroscopic evolution of 
SN 1997D for over 1 year after discovery. 

The photospheric spectra are dominated by a red con- 
tinuum and P-Cygni profiles of HI, Ball, Call, Nal, and 
Sell (see Fig. 4 of Paperl) . Spectral modeling showed that 
the unusual strength of the Ball lines compared with typical 
SNII is due to the low temperature rather than to overabun- 
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Dut = Dutch 0.9m+CCD Cam.; CT91 = CTIO 0.91m+CCD 
Cam. 

3.6 = ESO 3.6m+EFOSCl; 2.2 
NTT = ESO NTT+SUSI; Dan : 
EF2 = ESO 3.6m+EFOSC2 
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Table 3. Spectroscopic observations of SN 1997D 



Date 



phase* 
(days) 



inst. 



range 



(A) (A) 



14/1/97 


+33 


1.5 


3700-7570 


7 


15/1/97 


+34 


1.5 


3700-7570 


7 


16/1/97 


+35 


1.5 


3700-7570 


7 
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1.5 
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7 
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20 
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17 
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3.6 
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17+7 
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2.2 
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10 


2/2/98 


+417 


EF2 


3860-8030 


13 



relative to the estimated epoch of maximum, JD=2450430 
* See coding of Table 1, plus 1.5 = ESO 1.5m + B&C; CT4IR 
= CTIO 4.0m + IR spect.; CT1.5 = CTIO 1.5m + R-C spect. 
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Figure 3. Spectral evolution of SN 1997D. Wavelength is in the observer rest frame. The ordinate refers to the first spectrum (+33d), 
the other spectra are shifted downwards by 1 X 10~ 15 , 1.7 X lCT 15 , 2 X 10~ 15 , 3 X 1CT 15 , 3.9 X 10~ 15 , 4.2 X 1CT 15 , 5.6 X 10 -15 and 
6.6 X 10 — 15 respectively. For clarity some spectra have been multiplied by factors given in parenthesis. 



dance. However, the most striking property of these spectra 
is the very low expansion velocity of the ejecta. The minima 
of the absorption lines give an expansion velocity of about 
1100-1200 kms" 1 for Ha and somewhat larger (1500-1800 
kms -1 ) for Ball. Spect ral modeling suggests ph otospheric 
velocities between 900 ( |Chugai fc Utrobin 2000| ) and 970 
kms _1 (PaperI) and kinetic energies between 1 — 4 X 10 50 
erg, depending on the ejected mass. The spectral energy 
distribution in the optical-IR about 50 days past discovery 
is shown in Fig. ^. As mentioned above, the distribution is 
peaked in the optical window and very little flux is emitted 
at wavelength smaller than 5000 A. 

Two months later the Ha emission starts to dominate 
over the Ball 6497 line and becomes the more intense spec- 
tral line along with NaT D and Call IR lines. Starting 283 



days past maximum the spectrum is that of a typical SNII 
in the nebular phase but for the very low kinetic energy. 



2.3 Nebular line identifications 

The low expansion velocity of the ejecta of SN 1997D offers 
a unique opportunity to resolve and identify spectral lines. 
These identifications are presented in Fig. [sjfor the spectrum 
at +417d. 

The strongest feature after the prominent Ha is the 
7300 A doublet. Because of the close wavelength coincidence 
we associate this emission with the [Call] lines 7291-7323 A 
and definitely exclude the contribution of [Oil] 7320-7330 A 
doublet. We remind that the latter was instead identified in 



the late (about 15yr) spectra of SNe 1979C and 1980K (|Fc 
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are plotted for reference as a continuous, dashed and dotted line, 4000 5000 6000 7000 8000 

respectively. The IR spectrum is from Clocchiatti et al. (2000) rest wavelength (A) 



sen et al. 1999). Just blueward of this feature an unblended 
emission is identified with [Fell] 7155 A. Though sometimes 
questioned, this line was already identified in other SNII 
at late epochs (e.g. 1988A and 1988H, Turatto et al. 1993). 
Here for the first time we single out two other strong lines of 
the same [Fell] multiplet 14 (AA 7388 and 7453 A), confirm- 
ing this identification. Many other lines of [Fell] and Fell 
are identified in the spectrum. In particular, several lines of 
the Fell multiplet 42 and the [Fell] multiplet 19 are clearly 
visible. 

Another strong line is the emission due to [Mgl] 4571 
A. In contrast, H/3 is very faint if present at all, making the 
Balmer decrement > 45 which indicates collisional excita- 
tion. 

We could not find an unambiguous identification for a 
faint line measured at 6499 A (rest frame). For continuity 
with the photospheric spectra, we tentatively identify this 
line with Ball 6497 (mult. 2). To support this identification, 
we searched the spectrum for other Ball lines. It turns out 
that the strongest features (mult. 1 and 2) if present, would 
be blended with intense iron emissions, and therefore we 
cannot make any definite statement. 

More interesting is to analyze in detail the evolut ion o f 
the [OI] 6300-6364 A doublet. Spyromilio & Pinto (]l99l|) 



have shown how it is possible to derive an estimate of the 
OI density from the evolution of the intensity ratio of the 
two lines, which ranges between 6300/6364 = 0.95 when the 
nebula is optical thick to 3.06 in case the nebula is optically 
thin. In the three nebular spectra of SN 1997D, we measure 
6300/6364 = 0.9 ± 0.2 at +139d, 1.9 ± 0.1 at +283 d and 
2.65 ± 0.05 at +417d. Following Spyromilio & Pinto (jl99l|) 



we can estimate the oxygen number density at day +283 of 
1.5 ± 0.3 x 10 9 cm -3 . Noting that the FWHM of the lines 
(corrected for instrumental resolution) is 680 km s _1 and 
assuming uniform distribution, we estimate that the oxygen 
mass in the ejecta is ~ IMq. If the filling factor is less than 
unity (Spyromilio & Pinto 1991 estimate a filling factor of 
~ 0.1 for SN 1987A), the oxygen mass may be even lower. 



Figure 5. Line identifications on +417d spectrum. The wave- 
length is in the parent galaxy rest frame. 



Such an estimate seems more compatible with a low mass 
progenitor (see Chugai & Utrobin 2000) than with a high 
mass one, where the total mass of ejected oxygen is likely to 
be closer to 2Mq (Balberg et al. 2000) , but the latter cannot 
be ruled out in view of the large observational uncertainties 
involved in estimating the oxygen mass. Furthermore, the 
amount of oxygen left in the expanding envelope in the high 
mass progenitor case is also sensitive to the extent of fall- 
back on to the nascent black hole, so a low mass of ejected 
oxygen does not necessarily reflect the amount of oxygen in 
the progenitor. 



3 THE NATURE OF THE PROGENITOR AND 
THE COMPACT REMNANT 

In principle, the low luminosity of SN 1997D may be due 
to either a low 56 Ni mass or a low opacity to the 7-rays 
from radioactive decay. The fact that, for almost one year, 
the decline rate of the late time light curve of SN 1997D 
matches the decay rate of 56 Co is a strong indication that 
the envelope maintains thick to 7-rays and hence that the 
ejected 56 Ni mass is very low (0.002M Q ). To date this is 
the lowest measured value of the 56 Ni mass ejected in a 
supernova. The very low expansion velocity of the ejecta and 
the low explosion energy of a few 10 50 erg are confirmed by 
the spectrum taken 1.5yr after the burst. 

These observations single out SN 1997D as an extremely 
unusual type II supernova. As mentioned in the introduc- 
tion, there are two distinct types of progenitors that can pro- 
duce such an explosion - either a relatively high mass star (26 
M©; Paperl) or a relatively low mass one (8-10 Mq; Chugai 
& Utrobin 2000). The two models differ in their prediction 
regarding the nature of the compact remnant. For the mas- 
sive star model, significant fallback of matter is expected 
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because of the very low explosion energy and the deep grav- 
itational potential of such a massive star. In contrast, for the 
8-12 Mq model, the mantle of heavy elements and the sur- 
rounding He/H envelope are so extended that fallback is very 
small, even for explosions as weak as ~ 10 49 erg (Nomoto et 
al. 1982; Nomoto 1984, 1987; Nomoto & Hashimoto 1988). 
The significant amount of fallback in the high mass pro- 
genitor will most likely produce a black hole, while in the 
low mass scenario fallback is negligible and a neutron star 
is probably formed. As suggested by Zampieri et al. (1998a) 
and recently supported by Balberg et al. (2000), this dis- 
tinction is of special interest in the case of SN 1997D and 
in the high mass progenitor scenario, since it may offer the 
first opportunity of a direct observation of a newly formed 
black hole in a supernova remnant. 



3.1 Light Curve and Spectral Simulations 

We performed theoretical calculations of the late time light 
curve of SN 1997D for both high mass and low mass progen- 
itor models. The mass of the ejecta and the explosion energy 
of the high mass progenitor model (constructed from the 25 
Mq model of Umeda, Nomoto & Nakamura 2000) are M e j = 
18 M Q and E = 4 x 10 50 erg. The low mass model has M c j 
= 6 Mq and E — 1 x 10 50 erg, as adopted by Chugai & 
Utrobin (2000). All the calculations were performed with a 
spherically symmetric radiative transfer code, which solves 
a 7-ray transfer equation with one-energy group approxima- 



tion as well as a photon transfer equation ([wamoto et al 
200C| ). 

In both models, 56 Ni is distributed uniformly up to an 
expansion velocity of 800 kms -1 , as suggested by the spec- 
tral analysis (Paperl). We find that, in the high and low 
mass progenitor models, 99% and 94 % of the 7-rays are 
trapped at day ~ 400, respectively, and thus the optical 
light curves decline at a rate close to that of 56 Co, consis- 
tent with observations. At day ~ 400, the difference between 
the two models is too small to be detectable (the high mass 
model is only 5% more luminous than the low mass model). 
Significant deviations from the Co decay rate and differ- 
ences between the two models (~ 0.2 mag) are predicted to 
be seen only after day ~ 800. 

Models suggest that after the end of the plateau, which 
is the consequence of an inward-moving recombination wave 
in the H-envelope, the light curve of a SNII should be de- 
termined mostly by 7-ray deposition followed by immediate 
radiation of optical photons. Therefore the 56 Co rate should 
be a lower limit to the slope of the light curve. However, 
observations covering the period 60-200 days show that the 
decline rate was slightly slower than the 56 Co decay rate. 
This suggest that 7-ray deposition followed by instantaneous 
emission of optical photons might not be the only source of 
luminosity for SN 1997D at this epoch. This probably means 
that the photosphere has not disappeared yet. Indeed, spec- 
tra taken at 81d and 139d (Fig. 3) show the persistence of 
P-Cygni absorption profiles in many of the strongest lines 
(Nal D, Ball, Ha), even though the Call IR triplet is domi- 
nated by emission. Note also that the emission component in 
Ha is much stronger than the absorption component, while 
they are equivalent in Nal D. Although the contribution of 
net emission is significant, a thermal continuum might still 
influence the spectral shape. In addition we find that the 



B band, where there are no strong emission lines, declines 
much more slowly than the redder bands, all of which in- 
clude some strong emission lines. 

The presence of a photosphere at such advanced phases 
(its velocity, as measured from the P-Cygni profile of Ha, 
is ~ 1090 kms -1 ) might suggest the presence of an inner 
region with a steep density gradient, situated below the H- 
envelope. This fact might favor a high mass model and de- 
serves further quantitative study. The spectra taken at late 
times (283 days and 417days) show only emission lines and, 
in that period, the light curve decline rate becomes closer to 
the 56 Co decay. 

We tried to simulate the nebular spectrum of SN 1997D 
using a NLTE spectrum synthesis code that computes 7-ray 
deposition and heating, followed by cooling by forbidden line 
emission. The code is designed for Type la SNe, so H is not 
included. Fitting the observed Fell forbidden lines provides 
support to the finding that the Ni mass is of the order of 
0.002 Mq. Also, the velocity of the nebular lines is very 
small, about 800 kms -1 , confirming that the kinetic energy 
of the ejecta is low and that the small Ni mass is not mixed 
out significantly. 



3.2 Black Hole Remnant 

We now turn to discuss the implications of the late times 
measurements on identifying the nature of the compact 
remnant. The data indicate that the late time bolomet- 
ric light curve is consistent with the decay of 56 Co, and 
there is no evidence up to 400 days of a change in the de- 
cline rate caused by emission from the newly formed com- 
pact object. The bolometric luminosity at t = 417 days is 
7.8 x 10 38 erg s , sufficiently low that emission from a so- 
lar mass compact object accreting at the Eddington limit 
{L E dd ~ 2.6 x 10 38 A'/* erg s" 1 , where M* is the mass of 
the compact object in Mq and the accreting material is as- 
sumed to be hydrogen poor) should have already been de- 
tected through a deviation from a perfect 56 Co decay. Fur- 
thermore, the upper limit in the R band at ~ 700 days 
rules out emission at a luminosity level of ~ 10 38 erg s _1 
from the compact object at that time, and further strength- 
ens the conclusion that Eddington-limited accretion is not 
taking place. Therefore, quite independently of detailed the- 
oretical modeling, these observations in themselves tend to 
favor the BH formation scenario. 

The extremely low efficiency of converting spherical ac- 
cretion onto a black hole into radiative luminosity (for the 
accretion rates of interest) severely limits its ability to radi- 
ate at the observed rate. Hence, the absence of a signature 
from the compact object is qualitatively consistent with the 
presence of a black hole in SN 1997D. In fact, emergence of 
a black hole as early as t = 417 days would require a very 
high accretion rate over a relatively long time. The lumi- 
nosity from steady-state, hypercritical spherical accretion of 
helium-rich matter onto a b lack hole can be approximately 
expressed as (Blondin 1986): 



Lac 



1.44 xl(T 



Me 



Mq 



-1/6 



M 



Mq yr- 



5/6 



ergs 



"(1) 



so that to sustain a luminosity of a few xlO 38 erg s 1 at 
t = 417 days, the accretion rate would have to be a few 



© 0000 RAS, MNRAS 000, 000-000 



Supernova 1997D 7 



Mq yr _1 (depending on the exact mass of the black hole). 
Since the accretion rate is likely to be declining, it would 
have had to be at least as large prior to t = 417 days, im- 
plying that several Mq of material had been accreted in 
late-time fallback. For comparison, the best fit model for 
SN 1997D of Paperl has only 0.2 Mq of initially bound ma- 
terial available for late-time fallback (T. Young 1999, private 
communication). Indeed, the estimated accretion luminosity 
for 3 M Q black hole in SN 1997D at t ~ hundreds of days 
is a few times 10 36 erg s _1 (Zampieri et al. 1998a; Balberg 
et al. 2000), and emergence of the accretion luminosity is 
expected only at t~1000 days after the explosion. 

We note that if the bound material has sufficient angu- 
lar momentum, accretion should proceed via a disk, rather 
than through a spherical inflow. The disk would have to 
be very optically thick and advection dominated in order 
to accommodate the hypercritical accretion rates expected 
during the first years after the explosion (see Mineshige et 
al. 1997 for an hydrodynamic study of the properties of such 
disks up to ~ 1 day after the explosion). While the emis- 
sion properties of such advection-dominated, optically thick 
disks at hypercritical accretion rates are poorly understood 
at present, disk accretion onto a black hole can, in princi- 
ple, also be consistent with the current observations if the 
accretion efficiency is low. 



3.3 Neutron Star remnant 

The possibility that the remnant of SN 1997D is a neutron 
star can be examined in a similar fashion as the well studied 
case of SN 1987A. To date, the SN 1987A remnant has failed 
to provide any evidence for the presence of a neutron star 
formed in the explosion. 

First, we note that the luminosity at t — 417 days ap- 
pears to rule out the formation of a Crab-like pulsar in 
SN 1997D. The present magnetic-dipole emission estimated 
for the Crab pulsar is ~ 5 x 10 38 erg s _1 , and at an age 
of 417 days it was probably ~ 10 40 erg s" 1 (Kumagai et 
al. 1989; Woosley et al. 1989). If a sizable fraction of such 
emission is deposited in the supernova ejecta, it should have 
had an observable effect on the luminosity. 

If SN 1997D did produce a neutron star which is weakly 
magnetized or/and slowly rotating, this may have been ob- 
servable if it is accreting material through late-time fallback. 
Initially, the accretion rate is so high that such fallback is 
expected to be neutrino-cooled and to produce very little 
observable radiative luminosity, but this luminosity should 
gradually increase as the ex pansion of the e jecta causes a 
decline of the accretion rate (Chevalier 1989). 

In the case of SN 1987A, Houck & Chevalier (1991) 



found that the radiative luminosity should reach the Edding- 
ton limit about six months after the explosion, and remains 
at that value for several years. Following their approach we 
can make a similar estimate for accretion onto a putative 
neutron star in SN 1997D. At late times the accretion rate 
is expected to be ballistic (dust-like) with a secular time 
dependence of (Chevalier 1989; Colpi et al. 1996) 



M(t > t late ) 



1.2 x 10 



Po*o 



10 9 gcm- 3 s 3 



tlat 



-1/3 



10 days J 



-5/3 



Meyr' 1 . (2) 



Here po is the average initial density of the accreting ma- 
terial and to is the initial expansion time to = r/v (r is 
the radius and v the velocity of a gaseous shell), assum- 
ing initial homologous expansion at the onset of dust-like 
flow. ti a t e denotes the time at which the accretion flow set- 
tles into dust-like motion, and is roughly t; ate « 10 days 
for SN 1987A. Both the numerical coefficient and ti ate are 
weakly dependent on the details of the explosion. The quan- 
tity potg is related to the kinetic energy per unit mass £ and 
the amount of mass Mf, ounc j which is gravitationally bound 
to the neutron star at the onset of homologous expansion: 
p t 3 = 0.3 3/2 {3/4n)M bound £- 3/2 . 

The luminosity is expected to reach the Eddington limit 
roughly when the accretion shock approaches the trapping 
radius in the accretion flow. Thi s corresponds to an accre - 



tion rate of ~ 4 x 10~ 4 M Q /yr ([Houck fe Chevalier 199l|) 



Again, this value depends only weakly on the details of the 
explosion. Setting M ~ 4 x 10 _4 MQ/yr in equation (^j), we 
find that the Eddington-limited accretion phase is reached 
at time 



t, = 0.5 



3/5 



tlat 



10 days 



-1/5 



10 9 gcm _3 s 3 



3/5 



yr-(3) 



Unlike typical Type II supernovae which have a larger Co 
abundance, in SN 1997D this abundance is so low that the 
radioactive luminosity is comparable to the Eddington limit 
for a solar mass compact object as early as a few hundred 
days after the explosion. Hence, in SN 1997D, t, provides 
a reasonable estimate as to when the fallback luminosity 
will "emerge" above the emission of radioactive elements, 
causing the light curve to deviate from pure 56 Co decay. 

Assuming again initial homologous expansion and uni- 
form density throughout the post-explosion envelope, the 
amount of mass gravitationally bound to the neutron star 
can be roughly approximated as 

'M* \ / M ^ 1/3 



-Abound 



0.013 



poto 



10 9 gcm- 3 s s 



^10M Q 

2/3 

Mq 



U0 13 cm/ 



(4) 



where M and ro are the mass and radius of the envelope 
at the onset of expansion. Equation (^) shows that, even 
in extended progenitors (outer radius of several 10 13 cm), 
a small amount of material, of the order of 10 _2 Mq, can 
remain bound to the neutron star after the explosion and 
therefore be available for late time fallback. In fact, Mbound 
increases with increasing poto, i.e. with decreasing explosion 
(and kinetic) energy. This small amount of gas is dynami- 
cally negligible if compared to the mass of the neutron star. 
However, as shown by equation (^|), it is sufficient to estab- 
lish a considerable accretion rate at late times (~ 10 _4 Mq 
yr _1 at t ~ 1 yr) and to give rise to an Eddington-limited 
luminosity ~ 1 year after the explosion (eq. M). 

The post-explosion profile of the bound material in 



the 26 Mq progenitor model of Paperl has po^o 



1.5 x 



10 10 g cm -3 s 3 , for which the accretion rate would decline to 
4 x 10 _4 MQ/yr only at ~ 2.5 years after the explosion. Al- 
though this model involves forming a compact remnant with 
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a mass M J; 3Mq and therefore likely to be black hole, its 
parameters do allow for an accreting neutron star to remain 
hidden in its midst up to t — 900 days. On the other hand, 
the corresponding parameters for the low-mass progenitor of 
Chugai & Utrobin (2000), which is much more likely to have 
produced a neutron star, give po^o — 6.3 x 10 8 g cm -3 s 3 . 
This latter case is more similar to that of SN 1987A, and 
accretion onto the neutron star will have reached the critical 
rate 4 — 5 months after the explosion, and at t — 417 days 
would be radiating at the Eddington luminosity. The ab- 
sence of such luminosity in the light curve of SN 1997D at 
400 and 700 days could imply that the remnant is a black 
hole. 

Note, however, that the same reservations originally 
raised regarding SN 1987A apply here as well. First, ac- 
cretion luminosity would not be produced if the accretion 
flow is disrupted owing to dynamical instabilities. Fryer, Col- 
gate & Pinto (1999) have recently speculated that such an 
instability could arise for an accreting neutron star when 
the heavier elements begin to recombine, giving rise to a 
line-driven wind that eventually expels all the material that 
was initially bound. Furthermore, as noted by Mineshige, 
Nomoto, & Shigeyama (1993), if an accretion disk forms, 
the timescale for the accretion luminosity to decrease below 
the Eddington limit depends on the viscosity parameter, and 
further study is needed to investigate the evolution in this 
case. 



4 CONCLUSIONS AND DISCUSSION 

The analysis of the late-time photometric and spectroscopic 
data presented here highlight the peculiar character and 
unusual properties of SN 1997D. The very low expansion 
velocity of the ejecta (~ 1000 km s _1 ) and the very low 
abundance of 56 Ni (0.002 M ) make SN 1997D a unique 
supernova. The observations made until 400 days after the 
explosion (and an upper limit of the luminosity at 700 days) 
lend further support to the interpretation that this super- 
nova was the result of an unusually low energy explosion 
(~ 10 50 ergs). The available spectroscopic data show also 
that at t = 417 days the supernova ejecta have not yet 
started to interact with the CSM. 

The observed light curve of SN 1997D after day 300 
is consistent with pure radioactive heating, and there is no 
evidence of a significant deviation caused by other energy 
sources in the expanding supernova envelope. A neutron star 
formed from the collapse of a low mass progenitor would 
have produced emission at the Eddington rate if it were 
rapidly rotating and strongly magnetized (like the Crab), or 
if it were slowly rotating, weakly magnetized and accreting 
due to late time fallback. Such emission would have caused 
the light curve to deviate from pure radioactive heating de- 
cay by t = 417 days. Therefore, the presence of a neutron 
star remnant seems less likely. On the other hand, a spher- 
ically accreting black hole is indeed expected to be unde- 
tectable in comparison to radioactive heating during the pe- 
riod of observation, and hence it could not have revealed 
its presence (a black hole accreting hypercritically through 
a disk could also be undetectable if the accretion flow is 
advection-dominated) . The formation of a BH remnant is 
thus consistent with the available data. 



Several interesting implications can be drawn from the 
conclusion that a BH may be present in SN 1997D. First, 
the large amount of fallback needed to produce a remnant 
of ~ 3Mq would imply that the progenitor was a massive 
star. In this respect, it is interesting to investigate how the 
explosion energy of a massive star could be reduced to only 
a few 10 50 erg. We note that Iwamoto et al. (2000) proposed 
a scenario in which the core of SN 1997D had a small angu- 
lar momentum because the angular momentum of the core 
was transferred via magnetic-field interaction to the mas- 
sive H envelope. In this scenario the explosion of SN 1997D 
took place as a low-energy SNII instead of as a hypernova 
as in the cases of SN 1998bw and SN 1997ef, which had 
similar masses. Furthermore, if such low-energy supernova 
have a non-negligible rate (but are only under-detected), 
they could have a significant influence on the formation rate 
and mass distribution of compact objects, and perhaps even 
on nucleosynthesis abundances of heavier elements. 

Most notably, if the remnant is a black hole, SN 1997D 
may be the first case where obtaining direct observational 
evidence of the emergence of a black hole in the light curve is 
technologically possible. Such an observation would confirm 
the theoretical assessment that black holes are formed in 
supernovae (Timmes et al. 1996; Fryer 1999). As mentioned 
above, in the case of SN 1997D, the luminosity due to late- 
time spherical accretion onto the black hole may thus be ob- 
servable about three years after the explosion (Zampieri et 
al. 1998a; Balberg et al. 2000). The specific power-law time 
dependence of this luminosity should allow to distinguish 
it from any radioactive source and also from any possible 
emission caused by circumstellar interaction (which in any 
case has not been observed as late as 1.5 yr after the ex- 
plosion). We therefore re-emphasize the need for obtaining 
deep observations of SN 1997D in the very near future. 
We will focus our attention also on a more recent supernova, 
SN 1999eu, which shows many similarities with SN 1997D 
and for which we hope to be able to secure very late photo- 
metric and spectroscopic observations. 
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